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ABSTRACT 

We present multi-epoch high-resolution spectroscopy and photoelectric polarimetry of 
^ ' the long-period polar (AM Herculis star) QQ Vul. The blue emission lines show several 

distinct components, the sharpest of which can unequivocally be assigned to the illumi- 

■ nated hemisphere of the secondary star and used to trace its orbital motion. This nar- 
row emission line can be used in combination with Na i-absorption lines from the pho- 
tosphere of the companion to build a stable long-term ephemeris for the star: inferior 
conjunction of the companion occurs at HJD = 244 8446.4710(5)-h.Ex0^15452011(ll). 
The polarization curves arc dissimilar at different epochs, thus supporting the idea of 
fundamental changes of the accretion geometry, e.g. between one- and two-pole accre- 
tion modes. The linear polarization pulses display a random scatter by 0.2 phase units 
and are not suitable for the determination of the binary period. The polarization data 

Q i suggest that the magnetic (dipolar) axis has a co-latitude of 23°, an azimuth of —50°, 

( ' and an orbital inclination between 50° and 70°. 

, Doppler images of blue emission and red absorption lines show a clear separation 

■ between the illuminated and the non-illuminated hemisphere of the secondary star. 

The absorption lines on their own can be used to determine the mass ratio of the 
binary by Doppler tomography with an accuracy of 15% — 20%. The narrow emission 

' lines of different atomic species show remarkably different radial velocity amplitudes: 

^\ K = 85 — 130 km s^^. Emission lines from the most highly ionized species, Hen, 

originate closest to the inner Lagrangian point Li. We can discern two kinematic 
components within the accretion stream; one is associated with the ballistic part, the 
second with the magnetically threaded part of the stream. The location of the emission 
component associated with the ballistic accretion stream appears displaced between 
different epochs. Whether this displacement indicates a dislocation of the ballistic 
stream, e.g. by a magnetic drag, or emission from the magnetically threaded part of 
the stream with near-ballistic velocities remains unsolved. 

Key words: Accretion ~ cataclysmic variables - stars: QQ Vul- stars: imaging - 
polarization - Line: profiles. 



* Based in part on observations at the European Southern Ob- 
servatory La Silla (Chile) with the 2.2m telescope of the Max- 
Planck-Society 

t Visiting astronomer, German-Spanish Astronomical Center, 
Calar Alto, operated by the Max-Planck-Institut fiir Astronomic, 
Heidelberg, jointly with the Spanish National Commission for 
Astronomy. 



1 INTRODUCTION 

QQ Vul was discovered with the HEAO-1 low energy detec- 
tors in the 0.15 — 0.5 keV band and catalogued by Nugent 
(1983). Spectroscopic, photometric, and polarimetric obser- 
vations by Nousek et al. (1984) immediately following its de- 
tection led to the classification of this object as an AM Her 
type cataclysmic variable (or 'polar') with an orbital period 



© 1999 RAS 



2 A. D. Schwope et al. 



of 222.5 min. Nousek et al. (1984) found that the centroids 
of the prominent blue emission lines showed a huge systemic 
velocity, 7 ~ 500kms~^, when fitted with a circular orbit 
velocity equation, v = K sin 27r(0 — 0o) + 7. 

An extensive observational study using X-ray (EX- 
OS AT), UV (lUE), and optical spectroscopy and photom- 
etry was presented in Osborne et al. (1986) and Mukai et 
al. (1986). They found pronounced dips in the soft X-ray 
light curve caused by photoelectric absorption in the mag- 
netically confined stream (out of the orbital plane). In addi- 
tion, their spectroscopy revealed a systemic velocity 7 com- 
patible with zero, indicating substantial changes in the line 
emitting region with respect to the Nousek et al. (1984) re- 
sults. They also identified a narrow emission line component 
{peak 2) which they associated with the heated surface of the 
secondary star. 

In the same year, McCarthy et al. (1986) presented their 
analysis of moderate-resolution phase-resolved spectroscopy 
obtained in a state of reduced accretion. They identified a 
number of narrow emission lines and tentatively assigned 
one of them to the illuminated part of the secondary. Sur- 
prisingly, they found a stationary emission line component at 
redshift ~500kms~^ which remained largely unexplained. 

In 1986, Mukai & Charles detected the secondary star in 
near-infrared spectra of QQ Vul and, in 1987, the same au- 
thors were able to trace the Nal AA8183/94-doublet around 
the orbital cycle. Their data gave the first secure measure- 
ment of inferior conjunction of the companion star. 

EXOSAT observations of QQ Vul showed substantial 
changes of the X-ray light curve which were tentatively inter- 
preted by Osborne et al. (1987) as changes between a simple 
one-pole and a more complex two-pole geometry. This sce- 
nario has been supported by Beardmore et al. (1995) using 
more recent ROSAT X-ray observations. 

The lack of an ephemeris with sufficient accuracy ham- 
pers the unique interpretation of all the data collected in the 
past. In this paper, we will present an accurate spectroscopic 
ephemeris using both narrow emission lines from the illumi- 
nated hemisphere of the secondary star and photospheric 
absorption lines observed on several occasions between 1986 
and 1993. We show that both features are equally well suited 
for tracing the secondary. Our new spectroscopic ephemeris 
replaces older ones which were built using the recorded times 
of linear polarisation pulses or crosscorrelations of optical 
light curves. We will show that the linear polarisation pulses 
may occur at random phase (within a certain interval) and 
are not suitable, therefore, for the derivation of the binary 
period. 

We also present the results of polarimetric observations 
obtained between 1985 and 1988. These data are discussed 
in conjunction with the spectroscopic observations in order 
to establish a consistent picture of the accretion geometry 
of QQ Vul. Previous polarimetry was presented by Nousek 
et al. (1984) and by Cropper (1998), who derived the or- 
bital inclination to be 46° < i < 74° and i ~ 40° - 50°, 
respectively. 

This paper concentrates mainly on the analysis of the 
blue emission lines; in the companion paper {Mapping the 
secondary star in QQ Vut) by Catalan, Schwope and Smith 
(1999, henceforth referred to as Paper 1) we present a de- 
tailed analysis of the absorption and emission lines in the 



red spectral regime. Preliminary results were presented in 
Catalan et al. (1996) and in Schwope et al. (1998, 1999). 



2 OBSERVATIONS 

2.1 High-resolution spectroscopy 

QQ Vul was observed between 1986 and 1993 using the 2.2-m 
telescopes at La Silla and Calar Alto with Boiler & Chivens 
Cassegrain spectrographs, with the Calar Alto 3.5-m tele- 
scope equipped with the Cassegrain double-beam spectro- 
graph TWIN, and with the 4.2-m William Herschel telescope 
(WHT) at La Palma equipped with the double-beam spec- 
trograph ISIS. All the blue spectra covered the approximate 
wavelength range 4200-5000 A, where H-Balmer and Helium 
emission lines are prominent. The red spectra were centred 
on the Nal-doublet at AA8183/8194 A. FuU phase-coverage 
was achieved on all occasions. Some details of the spectro- 
scopic observations are given in Table |l[ All observations 
were performed under good weather conditions. Flatfield, 
bias, dark, and arc line spectra were taken frequently, and 
spectra of spectrophotometric standard stars were taken on 
each of the observing nights. We estimate the photometric 
accuracy of our spectra to be better than 30%. 

2.2 Photoelectric polarimetry 

Polarimetric observations of QQ Vul were carried out 
during several nights in 1985, 1987 and 1988 using the 
two-channel photometer/polarimeter of the DSAZ (Proetel 
1978) mounted on the 2.2m telescope at Calar Alto. In addi- 
tion, a rather short (2'') observation covering just one linear 
polarization pulse was carried out in November 1987 at the 
2.2m telescope at La Silla using the ESO-polarimeter PISCO 
(Stahl et al. 1986) . Most observations were done in unfiltered 
(white) light, hence the spectral bandpass was limited by the 
atmospheric cutoff and the sensitivity of the photomultipli- 
ers to ~3400 — 8500 A. An observation log is given in Table^ 
The sky was monitored regularly, allowing sky-subtraction 
from object measurements after fitting polynomials to the 
sky measurements. 

Using the DSAZ-polarimeter, the incoming signal is 
modulated by a rotating quarter-wave or half-wave plate, 
for the detection of circular (CP-mode) and linear polari- 
sation (LP- mode), respectively. In the CP-mode a linearly 
polarised signal can also be detected but with reduced effi- 
ciency. 



3 ANALYSIS 

3.1 Spectroscopic ephemeris of the secondary star 

The orbital mean spectrum of our 1993 observations (blue 
spectral range) is shown in Fig. |l|. The mean spectra ob- 
tained on the other occasions look very similar to that shown 
in the figure, although at a somewhat different brightness 
level. QQ Vul displays the typical features of an AM Herculis 
star in a high accretion state: a strong line of ionized Helium 
He II A4686, an inverted Balmer decrement, asymmetric vari- 
able lines with line wings extending to -^2000 km s~^, and 
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Table 1. Log of spectroscopic observations of QQ Vul 



Date 


Tel." 


Instr.*" 


Res. 


AA 


Tint 


No. spec. 


[Y/M/D] 






[A] FWHM 


[A] 


[sec] 




86/10/14-20 


LS22 


B&C+RCA 


3 


4180-5050 


600 


31 


86/10/18 


LS22 


B&C-fRCA 


1.5 


4540-4960 


625 


10 


88/06/10-12 


CA22 


B&C+RCA 


3 


4150-5050 


720/500 


22 


91/07/08-10 


CA35 


CTS-I-RCA2 


1.7 


4400 - 4940 


660 - 720 


38 


91/07/08-10 


CA35 


CTS-I-GEC 


2.2 


7550 - 8800 


660 - 720 


36 


93/08/23-26 


WHT42 


ISIS-fTekl 


1.5 


4230 - 5010 


300 


64 


93/08/23-26 


WHT42 


ISIS+EEV3 


0.7 


7940 - 8380 


300 


64 



" LS: La Silla, CA: Calar Alto, WHT: William Herschel Telescope La Palma; numbers follow- 
ing the observatory code denote the aperture diameter in units of 10 cm 

instrument codes: B&C - Boiler &; Chivens spectrograph, CTS — Cassegrain Twin Spectro- 
graph, ISIS — Intermediate-Dispersion Spectroscopic and Imaging System, the abbreviations 
following the instrument code denote the CCD chip used 



Table 2. Photoelectric polarimetry of QQ Vul 



Date Tel."" Instr.'' Mode'^ Filter Ti„t Ttot 

[Y/M/D] [sec] [hours] 



85/06/11 


CA22 


ZPP 


CP 


WL/RG630 


60 


3.9 


85/06/12 


CA22 


ZPP 


LP 


WL 


72 


4.1 


87/06/25-27 


CA22 


ZPP 


CP 


WL 


100 


8.2 


87/09/04 


LS22 


PISCO 


LP 


WL 


130 


2.0 


88/06/14+16 


CA22 


ZPP 


LP 


WL 


60 


5.0 


88/06/15 


CA22 


ZPP 


CP 


WL 


60 


4.2 



for coding of telescopes see Tab. |lj 

ZPP: two-channel photometer /polarimeter, PISCO: Polarimctor for Instru- 
mental and Sky polarization Compensation 
CP: circular polarimetry, LP: linear polarimetry 
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Figure 1. Orbital mean blue spectrum of QQ Vul in August 1993 



a pronounced Bowen blend of Clll/Nlll-lines between 4635 
and 4650 A. 

Individual emission lines consist of several compo- 
nents with different radial velocity amplitudes and different 
brightness variations. To give an example we show in Figj^ 
a grey-scale representation of the phase-folded, continuum- 
subtracted trailed spectrogram of the He ll A4686-line (1991 
observations). The ephemeris used for phase-binning is that 
given in Eq. 1. The trailed blue emission line spectrogram is 
shown together with a trailed, continuum-subtracted spec- 



trogram of the Na I absorption line doublet recorded simul- 
taneously. These red spectra are fully analysed in the ac- 
companying Paper 1 and are shown here to facilitate the 
interpretation of the blue spectra. 

The profile of the He-line generally shows great simi- 
larity to those displayed in HU Aqr in its high accretion 
state (Schwope et al. 1997). The most pronounced feature is 
a narrow emission line (NEL) with rather low velocity am- 
plitude. It moves parallel to the Na absorption lines, while 
its brightness is anticorrelated with those of the absorption 
lines. The latter are bright at inferior conjunction (blue to 
red zero crossing of the lines) , the NEL is bright at superior 
conjunction. The NEL clearly has a much lower radial veloc- 
ity amplitude than the absorption lines. Both the different 
brightness variation and the different radial velocity ampli- 
tude suggest that the two spectral features originate mainly 
from opposite hemispheres of the secondary star. This is 
demonstrated clearly by a combined Doppler map of both 
lines shown in Fig. ^. The underlying broader features in 
the Hell emission line complex with higher radial velocity 
amplitudes must then have their origin somewhere i n the 
accretion stream. They will be investigated in section 3.3.4. 
Here we concentrate on the narrow emission line which ob- 
viously originates somewhere on the donor star. 

We cannot confirm the finding by McCarthy et 
al. (1986) of a stationary narrow emission component at a 
velocity of '^500kms~^. 
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VELOCITY (km/sec) VELOCITY (km/sec) 

Figure 2. Trailed spectrograms of the simultaneously recorded Hell A4686 emission and Nal AA8183/94 absorption lines of QQ Vul (1991 
observations). The original spectra were continuum-subtracted and phase-folded according to the final ephemeris of Eq. 1. All data are 
shown twice for clarity 
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Figure 3. Doppler image (MEM-reconstruction) of the Hell 
emission and the Nal absorption lines recorded simultaneously 
in 1991. The Doppler image of Nal is shown as contourplot with 
iso-intensity lines at 20, 40, 60, 80 and 98% of the maximum inten- 
sity. The Doppler map of He II is shown as a grey-scale image with 
dynamic range between and 90% of maximum intensity. The 
overlay gives the size of a Roche lobe for a mass ratio Q = 1.75, 
the X marks the adopted centre of mass of the secondary star. 

The He ll A4686 emission lines in all our spectra ob- 
tained during different years were fitted by two or three 
Gaussians in order to separate the NEL from the stream 



Table 3. Fit parameters for sine fits to the Hell narrow emission 
and the Nal absorption lines 



Year 


Cycle 


7 


K 


to 


O-C 






kms^^ 


km s"-*- 


HJD" 


days'" 


Hell 


narrow emission line 








1986 


-11183 


-8±11 


118±11 


6718.4720(20) 


-7 


1988 


-7266 


-66±12 


144±13 


7323.7330(30) 


-H51 


1991 





-13±5 


116±6 


8446.4767(31) 


-H57 


1993 


5035 


2±8 


111±5 


9224.4817(19) 


-H19 


Nal absorption 


line 








1985 


-14396 


17±23 


209±28 


6222.0000(30) 


-f5 


1991 





-3±6 


271±8 


8446.47042(62) 


-6 


1993 


5028 


18±7 


228±10 


9223.3982(11) 


-HI 



"+244 0000 
* X 10-4 



emission. This was unequivocally possible only in the phase 
interval 0.25 — 0.60, at other phases the NEL is too faint 
to be clearly detected or strongly blended with another rel- 
atively sharp underlying component (around phases ~0.2 
and ~0.7). The centroid positions of the NEL were then 
fitted assuming a circular orbit velocity equation v = 'y -\- 
Ksm.2n{t — to)/ P, with fixed period P. 

Similarly, the Na I absorption lines were fitted by a dou- 
ble Gaussian with a fixed peak separation and the same line 
widths. Again the resulting radial velocities were fitted by a 
circular orbit velocity equation and the values given in Pa- 
per 1 are quoted here. Ellipsoidal fits gave smaller residuals 
than the straight sine fits (Paper 1); the influence on the 
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white light. All data were phase-folded and averaged using 
the ephemeris derived in Sect. 



3.1 



-20000.0 -10000.0 0.0 

ORBITAL CYCLE 

Figure 4. Residuals of events around the binary orbit in QQ Vul 
with respect to the linear ephemeris of Eq. 1. Meaning of symbols: 
o — NEL blue-to-red zero crossing; filled square - Nal blue-to-red 
zero crossing; X — Mukai et al. peak '2' component; ■ - linear 
polarisation pulses; * - X-ray absorption dips. 



times to of zero crossing, however, is negligible. We compile 
the results for all fit-parameters in Table H and include also 
those found in the literature (Mukai & Charles 1987, Nal 
measurement). 

Linear regressions for the times of blue-to-red zero 
crossing were performed for the emission and the absorp- 
tion features. They yielded identical results within the er- 
rors. The weighted linear regression for the combined Nal- 
and He 1 1- lines yields 

To(HJD) = 244 8446.47105 + X 0.15452011 (1) 

(48) (11) 

(the numbers in brackets are the uncertainties in the last 
digits). The epochs defined by integer numbers of E are re- 
garded as epochs of inferior conjunction of the secondary 
star in QQ Vul. All phases in both papers refer to the 
ephemeris of Eq. 1. The residuals of several features recorded 
around the orbital cycle, such as radial velocity zero cross- 
ings, centres of linear polarization pulses. X-ray absorption 
dips, are plotted in Fig. ^ with respect to Eq. 1. This dia- 
gram shows the good agreement of the narrow emission-line 
zero crossings with Nal absorption- line zero crossings. It 
shows furthermore that the 'peak 2' component of Mukai 
et al. (1987) which was suspected to originate from the sec- 
ondary star has nothing in common with the narrow emis- 
sion line from the companion investigated here. In addition 
it shows the variability of the phase of the X-ray absorption 
dip (the bars in the diagram denote the duration of that 
feature, not the phase error) and it shows finally, that the 
linear polarisation pulse appears at random phase between 
phase 0.33 and 0.50 and cannot be used in order to derive 
the orbital period. 



3.2 Analysis of polarimetric observations 

3.2.1 Overall characteristics of the light and polarization 
curves 

In Fig. P we show phase-averaged light and polarization 
curves with full phase coverage obtained in 1985 and 1988 in 



where phase zero indi- 
cates inferior conjunction of the secondary star. The data 
obtained in 1987 are very similar to those obtained one year 
later, but with poorer signal-to-noise. These data are there- 
fore not discussed separately. 

The 1985 data are quite similar to those published 
recently by Cropper (1998) and to those by Nousek et 
al. (1984) but show some striking differences from the 1988 
data with respect to the shape of the curves and the phasing 
of individual features. 

The light curve at all three epochs is double-humped 
with primary and secondary minima separated by a half 
orbital cycle and the primary minimum centred around 
^spoc — 0.9. In June 1985, maximum light occurred at the 
phase of the second hump whereas in all later data the first 
hump was stronger (see Fig. |l^). The phasing of the light 
curve does not show dramatic changes between 1985 and 
1988, quite unlike the circular and linear polarisation curves 
which do. 

Circular polarization is almost always negative with 
short excursions at low level towards positive values some- 
where in the phase interval (/)spcc = 0.30 — 0.55, and shows 
a pronounced depression at phase (jispcc = 0.80 and ^spoc = 
0.95 in 1985 and 1988, respectively. 

Linear polarization is always low, i.e. below 2-3%. 
The curves show one distinct linear polarization maximum 
(pulse) per orbital cycle. It is centred at phase (/)spcc — 0.42 
in 1985, and marks the end of the phase interval of van- 
ishing circular polarization in that year. In 1988, the pulse 
occurs at 0apec — 0.35, and marks the start of the interval 
of positive (or vanishing) circular polarization. Not so easily 
recognizable in the phase-averaged representation of Fig. H 
is a second fainter linear polarization pulse at (jispcc = 0.57 
observed in 1988 which coincides with a second circular po- 
larization zero crossing. The linear polarization angle dis- 
plays a large scatter in 1985, whereas the variation is much 
simpler in 1988. Meggitt & Wickramasinghe (1982) derived 
a simple relation for the derivative of the polarization an- 
gle at the time of the polarization pulse for the case of a 
point-like emission region, i/ip = cosi. The observed data do 
not allow us to determine a unique value for ip-p in 1985, 
contrary to 1988 where i ~ 60°. The fact that the value of 
■i/ip changes during the pulse in the 1985 observation means 
that the emission region must have had a more complicated 
shape on that occasion compared with the 1988 observation. 

The optical light and polarization curves of AM Herculis 
binaries are modulated by cyclotron beaming, eclipses by the 
secondary star and the accretion stream, and self-eclipses of 
the accretion region by the white dwarf itself. X-ray obser- 
vations on different occasions using the EXOSAT satellite 
have shown that there is no eclipse by the secondary star 
(Osborne et al. 1987). The polarimetric data in 1985 were 
taken simultaneously with an EXOSAT observation. The 
combined data show clearly that the depression in the circu- 
lar polarization at phase ^apec — 0.8 is due to an obscuration 
of the accretion spot by the intervening accretion stream. It 
is that part of the stream which is magnetically coupled and 
lifted out of the orbital plane. As has been pointed out by 
Cropper (1998), this requires the orbital inclination i to be 
larger than the colatitude of the accretion spot &a (measured 
with respect to the rotation axis). 
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Figure 5. Phase-averaged light and polarization curves of QQ Vul obtained in 1985 and 1988 in white light. Showrn are from top to 
bottom the brightness (in arbitrary units), the degree of circular and linear polarization and the polarization angle. The data are shown 
twice for clarity. Phase zero is the inferior conjunction of the secondary star (ephemeris given in Eq. 1). 



Around the phase of the X-ray and circular polarization 
dip the observer looks almost directly onto the accretion 
region (minimum angle between the line of sight and the 
surface normal at the accretion spot). At this viewing angle 
the optical flux and the circular polarization are reduced due 
to strong cyclotron beaming. The primary optical minimum 
(and circular polarization dip) is, therefore, caused by the 
combination of depolarization/absorption in the stream and 
cyclotron beaming. 

The secondary optical minimum at phase 0spcc — 0.4 
marks the phase when the accretion spot is at the limb of 
the white dwarf. Its phasing coincides with a sign reversal 
of circular polarisation and the occurrence of linear polari- 
sation pulses. On the basis of those clues, we interpret the 
pronounced secondary minimum in the hght curve seen in 
1988, which is even more pronounced in the 1991 and 1993 
data (see Fig. |l^, as caused by a partial self-eclipse of an 
extended accretion region. The EXOSAT X-ray hght curve 
of September 14/15 1985 shows a broad primary minimum 



at that phase which is being interpreted by us as the com- 
bined effect of the partial self-eclipse and foreshortening of 
an extended accretion spot. 

3.2.2 Modelling the 1988 polarimetry 

To determine the accretion geometry of QQ Vul we take 
the simple view of a pointlike accretion region in a dipolar 
geometry. A simplified sketch of the geometry is shown in 
Fig. 1^ Three vectors play an important role, the magnetic 
axis /I, the surface normal in the accretion spot s, and the 
field in the spot /. All vectors are shown projected onto an 
x' ~z plane centred on the white dwarf with z being parallel 
to the rotation axis. For an analogous representation of the 
geometry see e.g. Cropper (1989). The stellar co-latitudes of 
the dipolar axis, the spot and the field therein are measured 
with respect to the rotation axis and designated & with the 
corresponding subscripts ^, s, f. Similarly, the azimuths, x, 
of these vectors are measured with respect to the line join- 
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Figure 6. Variation of the linear polarization position angle in 
the 1988 data shown together with some models. The three model 
curves correspond to i = 70° (short-dashed line), 60° (solid line), 
and 50° (long-dashed line) and length of self-eclipse of 0.21 phase 
units. 



ing both stars. Conventionally the accretion stream in our 
geometry is deflected towards negative azimuth. Below we 
will discuss also the azimuth of the threading region, Xth 
(Fig. which is defined in the same reference frame. 

Previous determinations of the accretion geometry 
based on polarimetry were published by Nousek et al. (1984) 
and Cropper (1998). Nousek et al. used a point-like accretion 
spot and derived a stellar latitude of the accreting magnetic 
pole in the range 27° — 10°, and an orbital inclination of 
i = 46° — 74°. Cropper (1998) investigated extended ac- 
cretion arcs with constant temperature, magnetic field and 
density, and derived i ~ 40° — 50° and a co-latitude of the 



magnetic dipolar axis 



10°. 



The markedly different yet simpler behaviour seen in 
the 1988 data relative to all previous observations, called for 
a redetermination of the accretion geometry. Furthermore, 
with our extended data set, it is possible to find a solution 
that can encompass the results from both polarimetry and 
Doppler tomography. Previous analyses also lacked an ab- 
solute reference in the binary frame. We do not attempt to 
perform a least-squares fit of the polarimetric data, instead 
we try to globally optimize the model (including polarime- 
try, X-ray photometry and Doppler tomography). Cropper 
(1998), although using a more sophisticated model which ac- 
counts for an extended accretion arc, found a full inversion of 
the data (as in Potter, Hakala & Cropper 1987) objectively 
inefficacious and also only globally optimized his model. 

The following observed features are used for an esti- 
mate of the geometry: the length of the self-eclipse A(/>scif, 
the variation of the polarization angle, and the phase of the 
circular polarization dip. Variation of the dipolar axis /I, 
the inclination i, and the coupling radius then predicts a 
different run of the polarization angle. By integrating the 
equation of motion for one particle with given mass ratio 
Q, one can also predict its velocity pattern in the {vx,Vy)- 
plane (a schematic Doppler map) when projected with the 
assumed orbital inclination. The resulting image ca n then 
be compared with the observed Doppler maps (Sect. 3.3.4) 
thus constraining the azimuth Sth of the threading region. A 
solution is accepted if the predicted tomogram and the pre- 




Figure 7. Sketch of the dipolar geometry at the white dwarf 
surface. Shown are a cut through the white dwarf and a single 
dipolar field line projected on an appropriately rotated x' — z 
plane. Also shown are vectors indicating the orientation of the 
dipolar axis fl, of the surface normal in the accretion spot s, and 
of the magnetic field line in the accretion spot /, the latter being 
tangent to the field line at the footpoint. In a dipolar geometry s 
and / are always different. 



dieted curve of polarization angle are globally optimized and 
the phase of the X-ray/circuIar polarization dip is reflected 
well. 

The separation between the linear polarization pulses is 
0.24 phase units. The sign reversals of the circular polariza- 
tion are separated by 0.18 phase units. For our estimate we 
use the mean of both values, Ac/>aeif = 0.21. Both features 
are, in principle, measures of the same event, namely when 
the line of sight is perpendicular to /. The length of the 
self-eclipse of the accretion spot is shorter than this phase 
interval. The inclination is approximately 60°, as suggested 
by the slope i/)p of the polarization angle during the pulse 
(Fig. The values of i and Acpaeii constrain the possible 
values of the colatitude 5{ of the field line to a certain range 
(see Brainerd & Lamb, 1985, for technical details). The dip 
is centred at (/f>spcc = 0.95, corresponding to an azimuth of 
about -20°. 

Model curves for the polarization angle for the accept- 
able range of inclination angles, i = 50° — 70°, are shown 
in Fig. rJ together with the 1988 data. Deviations between 
model and observation in the phase interval 0.8 - 1.1 should 
not be taken seriously, since this is the phase of depolariza- 
tion in the stream and low intensity due to cyclotron beam- 
ing. The adopted values of the magnetic colatitude for the 
three model curves are S{ = 47°, 38° and 27° for i = 50°, 
60° and 70°, respectively. 

A globally optimised fit taking into account the different 
observational quantities was achieved with the following sys- 
tem parameters: mass ratio Q — M-1/M2 = 1.75, inclination 
i = 60°, orientation of the dipole axis 5^ = 23° and Xm = 
—50°, azimuth of the threading region xth = —20°. With 
these parameters the accretion spot emerges at &a = 33° and 
Xa = —40°, and the field line within the spot has 5f = 36° 
and Xf ~ —37°. In Fig. ^ a series of field lines around this 
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Figure 8. Model for the coupling of the accretion stream as seen 
by a hypothetical observer located at the white dwarf. The big 
circle represents the secondary star. The adopted orientation of 
the dipole /J is given by its azimuth Xm = —50° and its colatitude 
5u = 23°. 



particular one is shown as seen from the white dwarf. An 
observer viewing the binary at inclination i would be rotat- 
ing along a line of constant latitude I = 90° — i. Hence, the 
figure illustrates that the line of sight for a likely inclination 
of about 60° crosses the stream between x = —17°. . . — 32°, 



which corresponds to orbital phase 



= 0.91. . .0.95. At 



that phase the X-ray absorption dip was observed in the 
September 1985 data. The model thus reproduces the ob- 
served motion of the polarization angle, the observed phase 
of the X-ray /circular polarization dip, and it refl ects t he 
main features of the Doppler tomograms (see Sect. 3.3.4). 

Based on the analysis of the 1988 polarimetry we do 
not claim an accuracy of the orbital inclination i and the 
co-latitude (5/ better than 10°. One should bear in mind, 
however, that these two quantities cannot be varied inde- 
pendently, they communicate with each other via A(^s(,if- A 
higher inclination requires a lower co-latitude of the field 5/ 
and of the spot &a ■ We note in particular that our estimate of 
i based on polarimetry is in full agreement with the results 
of Paper 1, i = 65°, using high-resolution spectroscopy of 
the infrared Na lines. It is also in agreement with the analy- 
sis in Nousek et al. (1984), which is worth mentioning since 
their data look quite different. The somewhat lower values 
favoured by Cropper (1988), i — 40° — 50°, are likely to be 
caused by the large scatter of the polarization angle in his 
data which indicates a more complicated accretion geometry 
in 1983 - 1985. 

Our results imply certain values of the azimuth and co- 
latitude of the accretion spot, which are in good agreement 
with X-ray photometry (not used for the determination of 
the geometry so far). The value of Xs ~ —40° predicts a min- 
imum in the X-ray fiux at phase 0.39 due to projection ef- 
fects. Indeed, this is exactly the phase where the broad X-ray 



minimum occurred in the EXOSAT data taken in Septem- 
ber 1985 (Osborne et al. 1987). Although the X-ray and our 
polarimetric observations were not carried out simultane- 
ously, the simple shape of the X-ray light curve in 1985 and 
the simple polarimetric behaviour of QQ Vul seen in 1988 
are suggestive of a similar accretion geometry. On the other 
hand, our parameter combination of i and 6a predicts a com- 
plete self-eclipse of the point-like accretion region at that 
phase, i + 5s — 93° > 90°, which is not observed. However, 
this does not necessarily imply that the model is incorrect. 
For one, the uncertainties associated with our derived values 
are larger than 5°. More importantly, the accretion region 
is likely to be extended. 



3.2.3 Changes of the accretion geometry 

Up to this point, we have considered and modelled only the 
1988 observations. The shift of the X-ray absorption dip to- 
wards earlier phases in June 1985 and October 1983, when 
the X-ray light curve had a complex shape, suggests that 
threading of the stream occurs further downstream, i.e. far- 
ther away from the Li-point at an azimuth of Xth ~ —70°. 
The X-ray lightcurve of June 1985, on the other hand, shows 
a completely different shape from the September 1985 data, 
which cannot be explained by a simple migration of the 
accretion spot towards a different azimuth. The fact that 
phasing of the X-ray maximum in June 1985 occurs ex- 
actly at the phase of minimum flux in the light curve taken 
three months later suggests the presence of a second active 
accretion region on the opposite hemisphere (Osborne et 
al. 1987). The corresponding cyclotron emission from the 
second region must affect the polarization measurements 
and probably distorts the simple variation of the polariza- 
tion angle seen in 1988, which likely represents a one-pole 
accretion mode. This was already pointed out by Cropper 
(1998) who also detected a second polarization pulse at opti- 
cal minimum phase. We conflrm this likely presence of a sec- 
ond accretion region based on our observation of enhanced 
linear polarization at <^spcc = 0.0 in our 1985 data. 



3.3 Analysis of the emission lines 

3. 3. 1 Origin of the narrow emission lines 

Our initial analysis of the simultaneously recorded spectro- 
gram s of the He ll emission and the Na I absorption lines in 
Sect. 3.1 based on radial velocity fits suggests that the NEL 



and the absorption lines originate on opposite hemispheres 
of the secondary star. This impression is supported by a 
combined Doppler tomogram of both lines (Fig. ^ which 
was constructed using MEM-deconvolution of the trailed 
spectrograms. We had chosen the He ll emission line for this 
initial analysis because it is intrinsically the sharpest of the 
bright lines in the blue spectral regime. With a full width 
at haff maximum FWHM of 1.9 A (1991) and 2.0 A (1993) 
at phase 0.5, the NEL of Hen is just resolved. The NEL 
components of the H-Balmer lines at the same phase have 
a FWHM of 3.7 A, those of Hei have 2.5 - 4 A. The NEL 
of Mgll A4481 is comparable in width to that in HenA4686 
but has a much lower flux and thus reduced signal-to-noise. 

A comparison of the results from the fit to the radial 
velocity measurements and those from the Doppler map 
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Figure 9. Numerical experiments showing tlie effects of noise 
and sampling on the reconstructed Roche lobe. The input Roche 
lobe used for generation of synthetic spectra is shown in each 
panel, (a) Reconstruction without noise, infinite velocity reso- 
lution and velocity sampling 12kms~^, phase sampling 0.00556; 
(b) reconstruction with noise at observed level, velocity resolution 
80kms~^, phase sampling 0.0556; (c) reconstruction with noise 
at observed level, velocity resolution 180kms~^, phase sampling 
0.0556; (d) reconstruction of the observed (1991) spectra. In pan- 
els (b) - (d) the isocontours represent between 10% and 90% of 
maximum intensity in steps of 20%. The isocontours drawn in 
(a) are the 20% - 80% intensity levels in steps of 20% with an 
additional contour at 90%. 



(Sect. 3.1) for the same data set reveals an inconsistency be- 
tween the radial velocity amplitudes in the NEL. Maximum 
emission in the Doppler map occurs at t; = jnyj ~ 96kms^^ 
whereas the radial velocity fit gives v — 116 ± 6kms~^. 
We think that this inconsistency results from the reduced 
phase interval in which the NEL velocities can be derived 
with good accuracy. Thus, we will use radial velocity am- 
plitudes derived from the Doppler maps which are based on 
all the data instead of a restricted data set. To compute the 
Doppler maps we mainly used the code of Spruit (1998). 

The relatively sharp separation between emission and 
absorption features in the Doppler map shown in Fig. ^ re- 
quires us to address the following questions. (1) Is it possible 
to reconstruct the shape of the Roche lobe and thus deter- 
mine the mass ratio Q using Doppler maps of the Na I- lines 
alone? (2) Is it possible to measure the size of the Roche lobe 
of the secondary directly using the velocities of the illumi- 
nated and the non-illuminated hemispheres of the secondary 
star provided these could be unequivocally determined ob- 
servationally? Answering these questions, in particular lo- 
cating the Li-point, is essential for the next section's inves- 
tigation of the accretion stream using Doppler tomography. 
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Figure 10. Velocity ratio of the photocentres of the irradiated 
front- and the non-irradiated backside of the donor star due to 
our irradiation model compared with the observed velocity ratios 
of Nal absorption lines (271 km s"'^) and narrow emission lines 
(NEL) of different species (data from 1991 and 1993). The solid 
line is a spline fit through the filled dots which were computed for 
i = 90°, the open circles were computed for i = 60° (fit indicated 
by dashed line). 



3.3.2 An irradiation model for the secondary star 

To address both questions we computed synthetic line pro- 
files using an illumination model of the secondary star like 
the one used in MR Ser and HU Aqr (Schwope et al. 1993, 
Schwope et al. 1997). The procedure is described in detail 
in Beuermann & Thomas (1990). 

In brief, the irradiating source is located at Mi, it is 
assumed that the secondary star fills its Roche lobe, that no 
shielding of the secondary's surface by the accretion stream 
or an accretion curtain occurs and that a sharp separation of 
irradiated and non-irradiated parts of the secondary exists. 
We assume that both emission and absorption lines origi- 
nate from the photosphere, i.e. from the Roche lobe of the 
secondary. The shape of the Roche lobe is iterated for given 
Q. Line emission from a surface element on the irradiated 
secondary is taken to be proportional to the solid angle as 
seen from the source at M\. Radiation is assumed to be op- 
tically thick, hence foreshortening of surface elements has to 
be taken into account. For the non- illuminated parts of the 
secondary star we assume that an absorption line is formed 
whose strength is proportional to the size of the surface ele- 
ment and also take foreshortening of the individual element 
into account. Emission and absorption line profiles are syn- 
thesized using the full Roche geometry. Since the spread of 
the radial velocity over the secondary's surface is of the order 
of several hundred kms~^, we regard velocity broadening as 
the only broadening mechanism of spectral lines. In order to 
transform velocities from dimensionless binary coordinates 
into the observer's frame, we adopt the mass-radius relation 
by Neece (1984) for the secondary star. The often used em- 



© 1999 RAS, MNRAS 000, 0-|l| 



10 A. D. Schwope et al. 



pirical relation by Caillault & Patterson (1990) gives slightly 
different results (see below). 

As mentioned above, our model does not take into ac- 
count any kind of shielding of the secondary star by e.g. an 
accretion curtain. This was present at some occasions in 
AM Her and HU Aqr (Davey & Smith 1996, Schwope et 
al. 1997) and became obvious in these stars by a prominent 
left /right asymmetry of the Doppler maps. Such an obvious 
asymmetry is not present in the combined He/Na Doppler 
map of QQ Vul, the centres of light of both species lie at 
Vx = Okms~^ (Fig- ^)) i-e. indicating no or far less shielding 
in 1991 than in 1993. Shielding was present in 1993 (cf. Pa- 
per 1), with the effect that the radial velocity curve of Nal 
had a higher ellipticity, yet smaller amplitude. 

In order to answer question (1) we computed absorp- 
tion line profiles with the same spectral resolution, orbital 
phase sampling and noise level as our 1991 observations. 
In Fig. ^ we compare reconstructions of data with (a) infi- 
nite resolution and S/N, (b) resolution, sampling and noise 
as in 1991 (upwhm ~ 80kms~^ measured using arc lines), 
(c) sampling and noise as before but degraded resolution 
(upwHM ~ 180kms~^) and (d) with a reconstruction of the 
observed data. The Roche lobe used for generation of the 
data is shown in each panel. 

From that exercise we learned (panel a) that it is pos- 
sible to reconstruct the original shape of the non-irradiated 
hemisphere of the companion star using high-quality data 
(high resolution, low noise), not a big surprise. The recon- 
struction shown in panel (b) , which uses data of same quality 
as the observed data appears somewhat smeared. This map 
has already lost the arc shape of the inner contours of the 
reconstruction in panel (a). 

Due to the limited resolution of the data, the contours 
extend into the nominally dark, front-side hemisphere of 
the star. The whole map has a smaller FWHM than the 
reconstruction of the observed data (panel d). In order to 
match the FWHM of the map in panel (d) we had to de- 
grade the resolution in our simulation to values as high as 
Av = 180kms~^. The resulting map (c) is almost rota- 
tionally symmetric and no longer displays a half star as all 
other maps do. The high degree of similarity between maps 
(b) and (d) suggest that the Na I absorption almost entirely 
originates from the non-irradiated hemisphere of the com- 
panion star and that the size of the Roche lobe chosen for 
the simulations is a fairly good representation of the true 
size. Variation of Q increases the size and shifts the lobe 
along the u^-axis and makes the fit worse. 

Our guess is that the observed map is disturbed by 
some effect and is, therefore, not as sharply reconstructed as 
the simulated one. For example, emission from the nearby 
HellA8236 line and the corresponding curved continuum 
could easily contaminate the absorption lines. The range 
of Q, which gave acceptable fits (by eye) to the location 
and curvature of the contours of the observed Doppler map, 
is Q = 1.5 — 2.2 (for i between 50° and 71°) with a best 
agreement at about Q = 1.8. 

Figure ^ addresses question 2. We attempt to (fur- 
ther) constrain the mass ratio using the velocity ratio of the 
photocentres of the illuminated and non-illuminated hemi- 
spheres of the secondary star. Synthetic spectra computed 
for different mass ratio Q and orbital inclination i were used. 
The inclination only has a small effect on the velocity ratio. 



Theoretically, the velocity ratio Vh/vi is a steep function 
of Q. This - in principle - offers the opportunity to mea- 
sure the size of the Roche lobe by determining the velocities 
of the photocentres. We assume that, observationally, the 
velocity of the non-irradiated side of the secondary star is 
Vb = 271 km s"'^, the largest velocity amplitude of the Nal 
lines measured so far (1991 data). Above we have argued 
that the Doppler map is consistent with a complete deple- 
tion of the irradiated front side of Na absorption. Incomplete 
Na depletion would mean that the true Vb would be larger 
and the derived mass ratio be smaller. 

Using all NEL's radial velocity amplitudes with esti- 
mated uncertainties better than lOkms"'^ as potential trac- 
ers oivf, observed ratios Vb/vf were computed. These ratios, 
shown as dotted lines in Fig. |l^ are highly dispersed, from 
85kms~^ for Hel4471 in 1993 to 135kms~^ for Call in 1991, 
and do not allow us to constrain the mass ratio better than 
when using the Nal lines alone (see above and Paper 1). 

For Q ~ 1.85 and i = 65° (Paper 1), the predicted 
velocity Vf is 115 — 125kms~^, that of the Li-point is 
v-Li ~ 75 - 80kms"\ Vb/vf = 2.36 - 2.17. For these pa- 
rameters the H-Balmer lines are apparently the most suit- 
able tracers of w/. The Ca and Paschen lines are compatible 
with this solution within the accuracy of our velocity de- 
termination, while the Hen lines display too low velocities, 
i.e. they originate closer to the Li-point. Obviously these 
lines are formed in an extended quasi-chromosphere located 
above the photosphere. The He ionizing radiation is not able 
to reach the whole geometrically allowed hemisphere of the 
secondary. The data presented here are thus able to resolve 
(although marginally) these different layers. A quantitative 
understanding of the formation of the emission lines requires 
detailed modelling of non-isotropically X-ray irradiated pho- 
tospheres of late-type stars, a project within reach of current 
models (e.g. Hauschildt, Baron & AUard 1997), but not yet 
undertaken. 

3.3.3 Multi-epoch line flux variations 

For the analysis of emission lines originating from the 
stream, we make use of emission line light curves and 
Doppler maps, shown in Figs. and 0. The emission line 
light curves are based on the integrated line profiles of the 
trailed spectrograms shown in Fig. [13. We approximately re- 
moved the contribution of the NEL, which does not originate 
from the stream, by measuring its fiux at phase (jispcc = 0.5, 
adapting a simulated emission line light curve to that level 
and subtracting it from the total flux. These NEL-free light 
curves are also shown in Fig. For comparison, contin- 
uum light curves extracted from a line-free region between 
the two emission lines H/3 and HellA4686 are also shown. 

The continuum light curves are all similar in shape dis- 
playing a double-humped structure with the primary mini- 
mum centred on ^spcc — 0.95 in 1991 and 1993 and proba- 
bly somewhat earlier, (/)apec — 0.92 in 1986. As discussed in 
Sect. |3.2| this minimum is caused by cyclotron beaming and 
the small shift of minimum phase at the different epochs 
indicates a possible azimuthal shift of the emission region. 

The emission line light curves are also double humped 
with intensity maxima at phase 0.2 and 0.7. This behaviour 
is more pronounced in H/3 than in Hell and more obvious 
after subtraction of the NEL, which peaks at (jispcc = 0.5. 
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Figure 11. (a, left) Continuum and emission line light curves (Hell A4686, H/3) of QQ Vul in 1986 (solid lines), 1991 (histograms), and 
1993 (circles). The continuum flux is the wavelength-averaged flux over the line-free wavelength interval indicated in the top panel in 
units of 10~^^ ergcm"^ s~^ A~^. Line flux is the integrated line flux in units of 10~^* ergcm"^ s~^. (b, right) Emission line light curves 
after subtraction of the approximated contribution of reprocessed emission from the secondary star thus leaving emission only from the 
accretion stream. 



Pronounced minima occur in the phase interval 0.92 - 1.00, 
the minima in H/3 generally occur earlier by ~ 0.03 — 0.05 
phase units at all epochs. The minima in both lines occurred 
earlier in 1986 than in 1993. The light curves in 1986 are 
highly asymmetric with a pronounced first maximum which 
is more than double the height of the second one. All these 
observed properties are suggestive of optically thick radia- 
tion in the accretion stream, with the stream having a higher 
thickness in 1986 than at the two other epochs. 



3.3.4 Multi-epoch trailed spectrograms and Doppler maps 
of He 11X4686 

Also the trailed spectrograms and the Doppler maps are 
markedly variable between the three epochs. We show in 
Fig. |l^ the trails and maps of HeiiA4686. The data of the 
other emission lines look very similar as far as the velocity 
pattern (though not the brightness pattern) is concerned. 
Looking first at the trailed spectra, the narrow emission line 
(NEL) is the most obvious and the most pronounced feature 
in all three observations. The underlying emission from the 
stream is brighter close to the NEL-emission in 1991 and 
1993 than in 1986. Hence, during the two recent observa- 



tions the stream was brighter in the vicinity of the inner 
Lagrangian point Li at rather lower velocities than in 1986, 
when stream emission was bright at higher velocities. 

This is nicely illustrated too in the Doppler maps which 
are shown as contour plots below the trailed spectrograms in 
Fig. |l2[ Contour levels were chosen appropriately to empha- 
size stream emission (the inner contour corresponds to the 
maximum intensity level from the stream). In 1991 and 1993, 
there is some bright emission originating approximately at 
Li and extending down to velocities Vx ~ —500 km s"'^ at 
constant Vy. From there a tail stretches down to the lower 
left region of the maps which can be recognized at veloci- 
ties as high as {vx,Vy) = (-1200, -700) kms"\ Within the 
resolution of our data, the orientation of these high- velocity 
tails is the same in all maps. 

The initial ridge of emission between Li and {vx, Vy) ~ 
(—400, 50) kms~^ can be tentatively identified with emission 
originating from the ballistic stream, as seen as prominent 
features in HU Aqr and UZ For (Schwope et al. 1997, 1999). 
There are, however, some problems with this interpretation, 
since this ridge shows a displacement by Avy ~ 45kms~'^ 
between 1991 and 1993. If it were a purely ballistic stream, 
it would be fixed in the binary and always appear at the 
same Wy-velocity. In the following, we will use the term 'bal- 
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Figure 12. Trailed spectrograms (top) and MEM-reconstructed Doppler maps (bottom) of tlie Hen emission lines in 1986, 1991 and 
1993. All spectra were continuum-subtracted and phase-averaged according to the ephemeris given in Eq. 1. The dynamic range of the 
grey-scale representations is between and 150 X 10~^^ ergcm"-^ s~^ A"'^, the same in each panel. Isocontour lines drawn with solid lines 
represent emission at 10%, 20% . . . 90% of maximum emission originating from the accretion stream. For the 1991, and 1993 tomograms 
contours at the 5% and 15% level are shown in addition drawn with dashed lines. 



listic' stream for the bright ridge of emission connected to 
the secondary star seen in the 1991 and 1993 tomograms, 
although it might not appear at a ballistic velocity at all. 

We determined the location of the ballistic stream by 
fitting a Gaussian curve to the Doppler maps at ~ 
— 250kms^^. Maximum emission ocurred at Vy ~ 45kms~^ 
in 1991 and at Vy ~ Okms"'^ in 1993. For a first-order esti- 
mate we set the iiy-velocities at Li and at = — 250kms~^ 
to be the same. For our best guess of Q and i, the projected 
velocity Vy at Li is 75 — 80kms~^ (see above). Only for 
the lower limit value Q = 1.5, i — 70°, is the projected 
Vy — 51kms~^ as low as observed in 1991. Hence, on both 
occasions (1991 likely, 1993 definitely), the observed ballistic 
stream appears at smaller velocities Vy than predicted by a 
single-particle trajectory for the likely parameter combina- 
tion of i and Q. 



Emission which could be apparently associated to the 
ballistic stream was almost completely absent in 1986. At 
that epoch a bright spot of emission centred on {vx,Vy) — 
(—380, —100) kms~^ emerged instead. 

As a first step to gain a quantitative understanding of 
the location and extent of the different parts of the stream 
in the magnetosphere, we have computed models for the ac- 
cretion stream, representing it by a single-particle trajectory 
under gravitational and centrifugal influence. We accounted 
for the pull of the magnetic field on the partly ionized stream 
by a magnetic drag force decreasing exponentially as a func- 
tion of radius. This additional drag redirects the trajectory 
in a small region in physical space. It initially follows a bal- 
listic path and then follows a dipolar field line. The velocity 
component along the field is conserved in our computation 
when the particle latches onto the field line. Depending on 
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Figure 13. Model of the accretion stream in Doppler coordinates. 
The adopted orientation of the dipole ji is given by its azimuth 
Xjj, = —50° and its colatitude 5^ = 23°, the same as in Fig. ^ 
With this geometry the ballistic part of the stream is truncated 
at Vx — — SOOkms"-*^. The magnetically controlled part of the 
flow appears in the lower left quadrant. The straight line at the 
centre of the bundle of trajectories indicates the location of the 
observed stream. 



displacement towards smaller Wy-velocities of the 'ballistic' 
stream with respect to the single-particle trajectory. One 
possible explanation for this displacement is that a proper 
ballistic stream (not influenced by the magnetic field) ex- 
ists but that the observed emission is dominated by matter 
with slightly changed velocity (matter in the process of get- 
ting coupled) . Another possibility has recently been worked 
out by Sohl & Wynn (1999) who treated the accretion flow 
as a number of diamagnetic blobs which independently in- 
teract with the magnetosphere of the primary. This model 
predicts blob trajectories differing from single-particle tra- 
jectories and might be applicable to QQ Vul. 

The apparent absence of a ballistic stream in 1986 
and the occurrence of the bright spot at {vx,Vy) ~ 
(— 380, — 100) kms^^ remains unexplained so far. There 
might be two explanations for the missing ballistic stream: 
either matter couples on the field directly at Li or the stream 
was outshone by other emission components and therefore 
not recognizable at the phase and spectral resolution of our 
observations. We investigated the first possibility by tun- 
ing the magnetic field strength in our models to high val- 
ues (which mimics low accretion rates), so that no ballistic 
stream was formed. We then used different orientations of 
the dipolar axis in order to check whether the trajectory 
runs through the observed bright spot down to the lower 
left quadrant in the {vx,Vy) plane. Azimuth and latitude of 
the magnetic axis were varied over large ranges but no solu- 
tion was found. We conclude that the resolution of the 1986 
data and the brightness of the other components prevented 
us from detecting emission from the ballistic stream. 



the orientation of the magnetic field in the threading region, 
the re-direction of the stream may result in a rather large 
step in velocity space {vx,Vy) (see Schwope et al. 1997 for 
further details). 

A model which accounts for most of the observed fea- 
tures of the 1991 and 1993 tomograms of Hell is shown in 
Fig. ^ We simulated several different trajectories coupling 
at different radii by varying the magnetic field strength, 
which appears as a dimensionless tuning parameter in our 
calculations. The trajectories shown in Fig. |l^ in velocity 
space are the same as those shown in Fig. ^ in real space. The 
smaller the assumed value of the field strength, the longer 
the particle remains on the ballistic trajectory (coupling at 
larger azimuth) and the higher the velocity at the threading 
point becomes. Observationally, the ballistic part of the ac- 
cretion stream can be followed down to — — 500kms~^. 
It is perhaps a bit more elongated in 1991 than in 1993, 
but this is difficult to state exactly, because the Doppler 
map of 1991 has a poorer velocity resolution. A velocity of 
Vx — — 500kms~^ on the ballistic stream corresponds to an 
azimuth of the coupling region of about —25° . Coupling may 
- in principle - occur all along the ballistic stream between 
the Li point and that azimuth. The orientation of the axis 
of the magnetic field was chosen in accord with polarimetry 
(Sect. ^^). With the chosen field geometry the magnetically 
funnelled part of the trajectories stretches along a more or 
less straight line in the lower left quadrant of the velocity 
plane, almost exactly along the observed high-velocity ridge 
in the Doppler maps (indicated by a straight line in Fig. hsl) . 

This simple model cannot account for the observed 



3.3.5 A simple model for the emission line light curves 

Support for the view that the bulk of matter does not cou- 
ple at Li but somewhere downstream comes from some nu- 
merical experiments which were performed to reach a basic 
understanding of the emission line light curves. A single- 
particle trajectory between the L\ point, the coupling re- 
gion and the white dwarf was divided in 100 parts of equal 
length and each element was assigned an intensity value. 
We allowed for two different brightness values, depending 
on whether the irradiated or the non-irradiated side of the 
stream was seen. Emission from an element was assumed to 
be optically thick, hence the observed brightness of an ele- 
ment at a given phase was the intrinsic brightness times a 
foreshortening factor. The brightness b of the vertical stream 
(along the dipolar field line) was varied as a function of ra- 
dius with h (X r^'^ starting with some initial value in the 
threading region. This radial dependence accounts for the 
shrinkage of the effectively emitting area due to the con- 
verging field lines. Light curves were computed by coadding 
over all visible elements at a given phase taking the full 
Roche geometry into account so that eclipses of parts of the 
stream are recognized if the inclination is high enough. 

In Fig. we show the resulting model light curves for 
the ballistic part of the stream and the magnetically fun- 
nelled part. The central trajectory of those shown in Figs. ^ 
and ^ was used and a contrast of 1 : 0.7 between the irra- 
diated and non-irradiated sides of the stream was adopted. 
The light curves were computed for an assumed orbital in- 
clination of 65°. Both light curves were normalized to max- 
imum emission. Emission from the ballistic stream shows 



© 1999 RAS, MNRAS 000, 



14 A. D. Schwope et al. 




-0.5 0.0 0.5 1.0 1.5 



SPECTROSCOPIC PHASE 

Figure 14. Model of the stream omission as function of orbital 
phase for the baUistic part in (a) and the magnetically funnelled 
part in (b). 

emission maxima at 0spec — 0.3 and 0.8 which is simply an 
expression of the deflection angle of ^18° of the ballistic 
stream with respect to the line connecting both stars. The 
primary minimum at ^spec — 0.02 is caused by the combined 
effect of foreshortening and a partial eclipse of the stream 
by the secondary star. 

The model light curve from the magnetically funnelled 
part of the stream shown in Fig. [l4[ 3 is affected only by 
foreshortening. This causes the maxima of emission to ap- 
pear at spectroscopic phases 0.2 and 0.7 and the minima 
at phases 0.45 and 0.95. These phases are well in agree- 
ment with the observed maxima and minima of the emis- 
sion lines of Fig. These phasings and the absence of 
pronounced stream eclipses by the secondary star are sug- 
gestive of stream emission mainly originating from the mag- 
netic part of the stream (in 1991 and 1993) or the threading 
region (1986). The visual impression of the Doppler maps 
(1991, 1993) is in apparent contradiction to this statement, 
since they show emission associated with the ballistic stream 
as a prominent feature. We regard the contradiction as ap- 
parent only for two reasons. Firstly, as discussed above, at 
least part of the emission associated with the ballistic stream 
originates from matter which is coupled to the field and 
thus is expected to produce a light curve resembling that 
in Fig. Secondly, the surface brightness distribution of 
a Doppler map does not reflect the brightness distribution 
in real space. In particular, emission at high velocities is 
spread over a large region in a {vx,Vy) diagram. The summed 
flux of the high surface-brightness pixels in the region de- 
flned by {vx,Vy) = (—500, —40) km s~^ (lower left corner) 
and (120, 200) kms"^ (upper right corner), i.e. including the 
NEL from the secondary star, contributes only to about 40% 
(1991) and 50% (1993) of the total intensity in the maps. 
Hence, low surface-brightness emission from higher veloci- 
ties is dominating in the stream. 

The Doppler maps indicate that we can best differen- 
tiate between emission from the ballistic and the magnetic 



stream at phases 0.25 and 0.75, i.e. in the direction of the 
iia;-axis. Projection in these directions produces a relatively 
sharp spectral feature of emission associated with the bal- 
listic stream (although not as sharp as the NEL from the 
secondary) whereas emission from the magnetic stream will 
appear as a rather broad feature. Fitting the observed spec- 
tra (1991 and 1993) at these indicated phases with triple 
Gaussians (NEL plus medium width plus broad component) 
reveals that the medium width and broad base components 
carry almost the same flux at these phases on both occasions. 
If we then take into account, that even the medium width 
component is fed by emission from matter which has already 
latched onto the field it seems plausible that most of the ob- 
served line emission comes from the magnetically funnelled 
stream and/or matter which is in the process of becoming 
funnelled, and that emission from the ballistic stream plays a 
minor role in the formation of the emission line light curves. 
Hence, Doppler maps and emission line light curves tell us 
the same story. 



4 CONCLUSIONS AND OUTLOOK 

Using medium- and high-spectral resolution observations 
with full phase coverage of the long-period AM Herculis 
binary QQ Vul obtained on three occasions in 1986, 1991 
and 1993 and complementary polarimetry obtained between 
1985 and 1988 we studied the accretion geometry in the 
inner and outer magnetosphere. First of all, we derived 
a reliable spectroscopic ephemeris for inferior conjunction 
of the secondary star by combining Nal photospheric ab- 
sorption lines with narrow emission lines of reprocessed 
emission from the X-ray irradiated hemisphere, HJD(ro) 
= 244 8446.4710(5) x 0? 15452011(11). The arrival times 
of linear polarization pulses show a large scatter in an 
(O — C) diagram with respect to the new ephemeris and 
thus are disqualified as useful tracers of the orbital motion 
for the case of QQ Vul. 

Analysis of the linear and circular polarization curves 
favour an orientation of the (dipolar) magnetic axis of 
5f_t = 23° and Xm = —50°, and an orbital inclination between 
50° and 70°. The variable shape of the polarization curves, 
in particular the phasing of circular polarization dips, sug- 
gests that the azimuth of the coupling (threading) region 
undergoes significant migrations in azimuth. 

The suggested orientation of the magnetic dipole and 
of the coupling region is in good agreement with the re- 
sults inferred from Doppler tomography of bright emission 
lines, in particular those of Hen A4686. These show (in 1991 
and 1993) three distinct emission features: (1) the irradiated 
hemisphere of the secondary star, (2) emission which can be 
associated with the ballistic stream, and (3) emission from 
the magnetically funnelled stream. Emission, which appar- 
ently arises from the ballistic stream, appears at different 
places in the Doppler tomograms of 1991 and 1993. This 
suggests that it is formed by matter which is in the pro- 
cess of becoming threaded, in complete agreement with the 
emission line light curves, which show the major contribu- 
tion from the magnetic part of the stream. At one occasion 
(1986) the Doppler maps do not reveal any emission that 
could be associated with the ballistic stream. 

We explored whether it would be possible to use the nar- 



© 1999 RAS, MNRAS 000, 



Multi-epoch Doppler tomography and polarimetry of QQ Vul 15 



row emission lines from the secondary star in combination 
with photospheric absorption lines to constrain the mass 
ratio and the orbital inclination. This was not possible in 
QQ Vul because the radial velocity amplitudes of different 
atomic species were found to be widely different. Lines from 
high ionization species such as HenA4686 originate more 
closely to the Li than a reprocessing model predicts. 

The Doppler tomograms presented in this paper sug- 
gest that threading of the accretion stream (or parts of 
it) starts very soon or perhaps immediately after leaving 
the secondary star. The stream seems to be completely dis- 
rupted at an azimuth of about —30°. Future observations 
with higher spectral and temporal resolution of the Nal- 
lines (1) will allow the mass ratio to be determined with an 
accuracy of better than 10% by the straightforward applica- 
tion of Doppler tomography and (2) will reveal substructure 
in the Doppler maps due to e.g. star spots. 
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